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The [dual] string tension of a spatial 't Hooft loop in the deconfined phase of Yang-Mills 
theory can be formulated as the tension of an interface separating different Zn deconfined 
vacua. We review the 1-loop perturbative calculation of this interface tension in the contin- 
uum and extend it to the lattice. The lattice corrections are large. Taking these corrections 
into account, we compare Monte Carlo measurements of the dual string tension with per- 
turbation theory, for SU(2). Agreement is observed at the 2% level, down to temperatures 
O(10)T c . 

1 Introduction 



Because of asymptotic freedom, the running coupling g(T) in an SU(N) Yang-Mills theory 
becomes small at high temperature. However, a perturbative calculation of SU(N) thermo- 
dynamic properties faces two obstacles. g(T) only runs logarithmically, so that subleading 
orders of the perturbative expansion contribute significantly as the temperature is lowered 
towards T c , the confinement/deconfinement transition temperature. And infrared diver- 
gences prevent an analytic, perturbative treatment altogether beyond some order. Heroic 
efforts have been devoted to the calculation of the pressure to this maximum order g 6 pP . 
The expansion converges poorly. One may fear that this is a general feature, and that an 
accurate perturbative calculation of any thermodynamic property is only possible at astro- 
nomically high temperatures. The purpose of this paper is to show a counter-example to 
this pessimistic view: the spatial 't Hooft loop. 

The tension of a spatial 't Hooft loop has been calculated in perturbation theory to the 
first three non-trivial orders, and the expansion appears to converge fast [2 ED- We show here 
that this perturbative calculation agrees with non-perturbative Monte Carlo measurements 
in the SU(2) theory to a high precision, down to temperatures O(10)T c . 

The 't Hooft loop has not received the same attention as the Wilson loop in numerical 
simulations of SU(N) lattice gauge theories. There are several causes for this relative neglect. 
• First, as 't Hooft has shown 0], area law or perimeter law for Wilson and 't Hooft loops 
together is forbidden, in the absence of massless modes. As a result, in SU(N) Yang-Mills 
theory, an area law is only observed for the spatial 't Hooft loop (dual to the temporal 
Wilson loop) above the deconfinement temperature T c . The associated dual string tension a 
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Figure 1: How to create a 't Hooft loop <9£ in the (x,y) plane: (left) in each (z,t) plane 
intersecting E, multiply by a center element one plaquette with coordinates (zo,to); (middle) 
equivalently, choose (20, io) in the corner; (right) equivalently, multiply by a center element the 
link Ut at the boundary. Thus, a 't Hooft loop of maximal size (E intersects all (z, t) planes) is 
equivalent to twisted boundary conditions for the Polyakov loop. 



serves as an order parameter for deconfinement |SJ|B]. The deconfined phase has traditionally 
received less attention than the confined phase, although the situation has been changing 
with the experimental search for the quark-gluon plasma. 

• Second, the intuitive, "physical" meaning of a spatial 't Hooft loop has not been widely 
appreciated, and it is often considered as an "exotic" observable of marginal interest. The 
fact is that a 't Hooft loop enforces a Zm interface in the Euclidean system, and that the 
dual string tension a is nothing else but the interface tension (up to a conventional factor 
T). 

This is easiest to see on the lattice. Starting from an Euclidean lattice of size L x x L y x 
L Z X Lt, with lattice spacing a (Li = Nia) and periodic boundary conditions in all directions, 
let us construct a rectangular 't Hooft loop W(dT,) in the (a;, y) plane. The contour <9E, which 
we take rectangular for simplicity, is the boundary of a surface E, both on the dual lattice. 
Let us take for E the minimal, planar surface, with coordinates (zo,to). Each plaquette of 
E is dual to a (z,t) plaquette. These (z,t) plaquettes with fixed coordinates (zo,io) form a 
"stack" 'P(E). The 't Hooft loop expectation value is then 

~ - = fVUcM-PEu P eV(t)(l ^ReTr (Up) - PY.u^v^ ~ I^r U P )) 
{ ( )} ~ fVUeM-PEu P (l-7f^U P )) 

(1.1) 

where £ € Zn (for SU(2), it amounts to flipping the coupling j3 — * — ff). Another choice of 
surface E leaves (W(dtl)) invariant, as can be shown by a change of variables in (|1 . 1|) : only 
the contour 9E matters. The dual string tension is then defined by 

a= lim --log(W(dt)) (1.2) 

A(8£)->-oo A 

where A is the minimal area bounded by 9E. We can then consider a 't Hooft loop of 
maximum size L x x L y , equal to that of the lattice in the x and y directions. In that case, 
the stack 'P(E) contains one plaquette in every (z,t) plane. Following Fig. for each (z,t) 
plane, we can move the special plaquette to a corner, then absorb the center element £ into 
the boundary link Ut 1 . As a result, the time-like links satisfy Ut(x,y,z + L z ,t — I) — 
(Ut(x, y,z,t = 1). The same happens for the Polyakov loop P(x, y, z) = FJ t U t (x, y, z, t) (the 
product of time-like links at position (x,y, z)), as 

P(x,y,z + L z ) = (P(x,y,z) . (1.3) 

Therefore, a 't Hooft loop of maximal size in the (x, y) plane is equivalent to twisted boundary 



1 The center element £ could be absorbed in U z instead, leading to the same free energy for a system with 
twisted boundary conditions for the z-like "Polyakov" loop. 
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conditions for the Polyakov loop, and 



(W(L X , Ly)) 




(1.4) 



where the numerator and denominator are the partition functions of a system with ordi- 
nary action, but with boundary conditions respectively twisted (by £) and periodic in the 
z-direction for the Polyakov loop. Taking the logarithm on both sides, dividing by L x x L v , 
and taking the thermodynamic limit, one recovers, on the left-hand side, the dual string 
tension tr, and on the right-hand side, the interface free energy per unit 2-dimensional area, 
i.e. the reduced interface tension, conventionally written as o~/T in a rather confusing no- 
tation. The dual string tension is identical to the reduced interface tension. Therefore, all 
the old studies of SU(N) interface tensions, both numerical JJj and perturbative [B], can be 
re-labeled as 't Hooft loop studies. It also becomes intuitively clear why a vanishes below 
T c : the Polyakov loop becomes disordered as the Zn center symmetry is restored, and the 
interface tension, i.e. ct, vanishes. 

• Third, the numerical study of the 't Hooft loop has been considered more difficult than 
that of the Wilson loop. The reason is an "overlap" problem. In the 't Hooft loop ex- 
pectation value ea. l|1.4|l . the two partition functions Z t ^ c and Z p ^ c are physically different. 
Gauge configurations which dominate the integral in the numerator Z t \, c contain an interface; 
configurations which dominate in the denominator Z p t, c do not. These two ensembles have 
little overlap, and importance sampling with respect to the denominator fails. This tech- 
nical problem can be approached in various ways, all of which entail multiple simulations 
and increased computer cost. Decisive progress was achieved with the "snake" algorithm 
[Bill] which factorizes the ratio Ztb c /Z p b c into N x x N y factors, each of 0(1), which can each 
be estimated by an independent Monte Carlo simulation. In each factor, the area of the 
interface is increased by one elementary plaquette a 2 . Finally, it was realized in 10 that a 
single factor converges to exp(— aa 2 ) in the thermodynamic limit. This brings the computing 
cost of a on a par with that of the string tension a of the Wilson loop (actually, we will 
see that the force a between two dual charges can be computed with a constant accuracy 
independent of their separation, in sharp contrast with the ordinary string tension). Contact 
with perturbation theory for a(T) becomes feasible at high temperature, as we will show. 

One surprise comes on the way. We expect our lattice measurement of a(T) to be affected 
by discretization errors 0{a 2 ), as for any bosonic theory, i.e. 



where c\ starts with 0(g 2 ). Here, c\ starts with O(g ), which can be traced to the non- 
perturbative nature of a. In this paper, we calculate the complete correction (<?(T, a) /<r(T, a = 
0) — 1) for g = 0. It turns out to be large and quite different from 0(a 2 ) for the lattice sizes N t 
accessible to current numerical simulations. In fact, it is not even monotonic as a function of 
Nt- The expected behaviour ~ c/N 2 , c <C 1, is only recovered for N t >8. These large lattice 
corrections obscure the analysis of lattice data, and mar the continuum extrapolations of 
old SU(3) studies of the critical interface tension [7j. Our 1-loop calculation of these lattice 
correction factors, shown in Table I and Fig. 0] should be of general use. 

Finally, we are in a position to compare our lattice Monte Carlo results for SU(2) with 
continuum perturbation theory 0^212 Ej. Remarkable agreement is seen (cf. Fig. 0), at 
the 2% level, down to temperatures O(10)T c , with perturbation theory at 0(g 2 ). 

Our paper is organized as follows. In Sec. II, we recall for completeness the perturbative 
calculation of the interface tension of Ref. [5], then show how this calculation is modified 
on the lattice, and extract the lattice correction factors Table I. In Sec. Ill, we compare our 
Monte Carlo results with the perturbative ones. Conclusions follow. 



a(T, a) = <j(T, a 



0)(1 + Cl a 2 + ...) 



(1.5) 
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2 1-Loop Calculation of the Interface Tension 



2.1 Continuum Derivation 

Following [Sj, we calculate the Z2 interface tension, which is derived from the effective action 
of a kink interpolating between the two SU (2) vacua. The effective action is the sum of the 
classical action plus a quantum term, obtained by integrating out fluctuations at 1-loop order. 

The calculation is done in euclidean space-time at temperature T . The euclidean time xo runs 
from to 1/T. The system has spatial size L x x Ly x L z , which is taken to infinity in the end. 

A Z(2) interface along the z-direction is constructed by assigning one vacuum to z — 
and the other to z = L z , and minimizing the effective action subject to these boundary con- 
ditions. By definition, the interface tension a is equal to the action of this interface divided 
by the area of the plane L x L y in the thermodynamic limit (a is often called the "reduced" 
interface tension; one can equivalently consider the (full) interface tension op, which is the 
interface action divided by the transverse volume L X L V /T. Of course, a = ap/T). 

We consider a gauge field A which is non-trivial only in the xo-direction: 

2ttT 

K = G ^ > Co = — g, Ci = i= 1,2,3 (2- 1 ) 
where T3 is the diagonal generator of SU(2) (T3 = Then the Polyakov loop is 

P(x) = -Tr 
K 1 2 

The trivial vacuum is at A'r = q = with -P(x) = 1. The Z(2) transform of the trivial 
vacuum occurs for q = 1, with P(x) = exp(j7r) = — 1. We can therefore introduce an inter- 
face by choosing q in (|2.1|l to be a function of z and fixing q(0) = and q(L z ) — 1. Then 
P(x) = 1 at z = and P(x) = -1 at z = L. 



Pexp ig / dxoAn (x) 



= cos(7rq) 



(2.2) 



We need to justify why the path of minimal action between the two vacua can be brought to 
the form Q2.ip . with q = q(z). That is because a constant background field can be brought 
to the diagonal, Cartan subalgebra by a global gauge rotation. But for SU(2), there is only 
one diagonal generator. 



The classical action is 



S C \A) 



alx / d 3 x -Tr [G^] 



(2.3) 



with the field-strength tensor G^ v — d^A v — d^A^ 
classical action is 

/ 2 2 
S = —^TL x L v ir 

g 2 



L, 



iglA^Au], For the field in (|2~T|l the 

2 



dz 



— {z) 

dz 



(2.4) 



Minimizing the action subject to the boundary conditions q(0) — 0, q(L z ) = 1, yields 
q(z) = z/L z . There is no true interface since the action vanishes like 1/L Z when L — > 00 is 
taken. But this is to be expected: classically there is no energy barrier, since the action is 
constant for any value of q. 



This degeneracy of the action with respect to q is broken by quantum effects. To show 
this, we calculate the action in the presence of the (classical) background field defined in 
1)2. To do this, we first assume q to be constant in space-time. The neglect of gradient 
terms will be justified a posteriori. 
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C Y m + Cgf = Tr 



S\x-y) (2.7) 



With = A C J + A q ™ the lagrangian consists of the following Yang-Mills term plus gauge 
fixing term (background field gauge) and Faddeev-Popov ghost term 

£ — &ym + &gf + £fpg 
£ym = \Tr[Gl u } 

Cgf = \tt[F*] 

C fpg = Tr[nPr,] (2.5) 

where the gauge fixing term is F = D^AJ^. The operator P in the ghost part is defined as 
P = ~2g^f^-, with a gauge transformation A. Also we define the covariant derivative in the 
adjoint representation for the background field: 

D$ = d li -ig[A*, } (2.6) 

The whole lagrangian quadratic in the fields is: 

Al u {x) (-S^Dl + (1 - j)DfD<^j Ar(y) 
C FPG = Tr [fj(x) (-2D 2 cl ) V (y)] 5 4 (x - y) (2.8) 
Using Feynman gauge (£ = 1) it follows: 

£ = £yM + £-GF + £-FPG 

= ±A%(x)M%{x,y)A%(y) + fj a {x)r ab (x,v)rn,{y) (2.9) 
T ab {x,y) = 5\x-y) [-D 2 cl ] ab (2.10) 
M%(x,y) = 8, v 5\x-y) [-D 2 cl ] ab (2.11) 

To obtain the quantum effective action to one-loop Sf u we have to integrate out the quantum 
fields A qu ,fj,r] 

e~ s " U = J VA qu VfjV V e- s 

= J VA^VffDnexp -\J d 4 xd%All(x)M;i(x,y)Al u b (y) 
-J d 4 xd 4 yfj a (x)J : ab (x,y)r] h (y) 

^Sl u = \Tr [log (M)]-Tr [log (?)] (2.12) 

To calculate such traces it is helpful to diagonalize the 'matrices' M and T by passing to 
momentum space. We then need the eigenvalues of the matrix 

/ k 2 + fc 2 + {2nTq) 2 2i(2TrTq)k \ 

-D 2 cl (k)=\ -2i(2TrTq)k k 2 + k 2 + {2nTq) 2 (2.13) 
V k 2 + k 2 ) 

At finite temperature, ko — 2nTn with n integer, and, using the definition 

k ± = 2nT(n±q) (2.14) 

it follows 

Ai = k 2 + [2nT(n + q)] 2 = k 2 + k\ 

A 2 = k 2 + [27rT(n - q)] 2 = k 2 + fc 2 (2.15) 
A 3 = k 2 + fc 2 (independent of q) 
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Figure 2: Left: The 1-loop effective potential is the function shown in the interval [0,1], then 
periodically repeated. Right: the kink solution for the effective potential, as a function of 

This shift in momenta is a typical effect of a constant background field. 

The sum over k+ is equal to the sum over fc_, therefore the quantum action reads, after 
substitution of (|!HU|) and lEUP into (E~I2>: 

s f v = v E /S 2 lo g( k2 + fc +) ( 2 - 16 ) 

n— — oo ^ ' 

where terms independent of q (eigenvalue A3) have been dropped. 

Now we need to perform the sums. To do this, wc first differentiate with respect to q. 

aqgu TV +00 r T +00 

^- = 5 E ^ ^ = - 8TW E (» + «) I" + 9l ( 2 - 17 ) 

'M— — 00 "t" n— — 00 

The divergent sum is interpreted using zeta-function regularization, with £(s, q) = X)n^o( n + 
q)~ s . Noting that Y^t=-oo( n + ?) I n + ?l = C(~ 2 , g) - ({-2, 1 - g), and using the iden- 
tity C(s, g) = ' B ~ s a +^ g - ) for s < 0, where £? is the Bernouilli polynomial, one gets = 
^T 3 Vn 2 [g 2 (l — q) 2 ] ■ Integrating back with respect to q, one obtains: 

Sl u = ?-T 3 Vn 2 [q 2 (l-q) 2 ] (2.18) 

So far we treated q as if constant, but actually we want it to be a function q = q{z) of z. 
The reason why we can still make use of the action S\ u for constant q is, that A cl varies 
slowly with z, and at leading order this variation can be neglected (gradient expansion), 
as we will see in the final kink solution (|2.21|l . Now we reintroduce the z-dependence of 
q — > q(z) and replace the size of the system in the z-direction L z (V — L x L y L z ) by an 
integral L z — > f^* dz. 

ST = \t 3 L x L^ 2 [ * dz[q(z) 2 (l-q(z)) 2 ] (2.19) 
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Now the classical degeneracy in q has been lifted by quantum effects and the minima (vacua) 
are at q = and q = 1 (see Figure |2 left). Since the sum l|2.17|l is periodic in q and hence 



6 



invariant under shifts q — > q + k for any integer k, this is also true for the quantum action. 
Thus q in Sf u (|2.19|l is defined modulo one. 



From (|2.4(1 the classical action is 



S — —TL x L y TT z 



dz 



dq ' 

Tz {z) 



Now introduce the dimensionless coordinate z' — cz with c = »J fgT in both S and Sf . 
The complete effective kink action then becomes: 

PI 2 



S 



kink 



s cl + Sf 



3 5 



I 1 L X LyTT 



L' 



dz' 



dq 
~dz~' 



(2.20) 



Defining the Lagrangian L = {^r\ + V(q), V(q) = q 2 (l — q) 2 , the Euler-Lagrange equation 

of motion gives -jp = y/V{q). The solution which satisfies the boundary conditions q{— oo) = 
and g(+oo) = 1 is 



q(z = z'/c) = — 



(2.21) 



This solution is called a kink (see Figure|21 right). The width of the kink is c 1 ~ which 
is much larger than the scale h of our system since j« 1. Hence it was indeed justified to 
neglect gradient terms ^§ : they only contribute at the next order in g 2 . 



Thus, we finally have for the effective action of the kink between the two vacua 



4 " 2 T \ L 



3%/6 g 

Our interface tension a is the kink action per unit planar area, i.e. S kmk / L x L y : 



a(T) 



4tt 2 T 2 



sVe g 



(2.22) 



(2.23) 



This formula superficially differs from the corresponding one in [H] by a factor T , but this is 
due to our different definition of the interface tension (S kmk = sc t v uf in |U)- 

2.2 Corrections on the Lattice 

When measuring the interface tension on the lattice, one may expect that the measurements 
of (N t ) (N t number of temporal sites) approach the perturbative continuum form (|2.23l) 
for g — * 0. Somewhat surprisingly, this is not the case. Multiplicative correction factors 
are needed, which are large for all Nt accessible to numerical simulations. Only for very 
large Nt do these factors go to 1. In \T2\ , Nathan Weiss gave the main expression leading to 
these correction factors, but without giving details about the calculation. This is done here 
for completeness. Moreover, we show how well these results describe the measurements we 
obtain from Monte Carlo simulations. 

Our derivation on the lattice is very similar to the one in the continuum and in the end, we 
compare the obtained effective potentials from both calculations. Actually the only thing to 
do is to find a lattice expression for the inverse propagator D 2 ,, which took the form H2.13(l 
in the continuum. 
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A four-dimensional lattice of size aN s in the spatial directions and aN t in the temporal 
direction is used. 



We calculate the classical action first. The Wilson action is 

= /?£(!" ^ReTrC/p) 



(2.24) 



Only plaquettes in the zt-plane are different from 1. Such a plaquette at position z takes 
the value 



with the forward derivative A{ . The trace of the plaquette is 

1 / vr 



2 cos ( a— A(q(z) 
Placing this into the action we have 

1 / 7T 



2 1- 



Pzt 



2 m 



^=/?e 5( a F A ^)) -iit- 2 - 2 !: (^w) 



g 2 N t 



(2.25) 



(2.26) 



(2.27) 



where t means sum over plaquettes in the zi-plane and fi = p-. 

Mimicking the steps from the continuum derivation, we add quantum fluctuations to this 
action. The calculation on the lattice yields for the momentum space propagator (|2.13|) 
(k — *• — where k' is now a dimensionless quantity) 



( E„4sin 2 (4)+4sin 2 (^)cos(fc^ 



-D 



2L 



2 

-2i sin 



V 



in(^) sin(fc^) 




2zsin(ff) sin(fc^) 
E,4sin 2 (^)+4sm 2 (f)cos(^) 



E,4sin(^) 



The eigenvalues are 



s * 

/ 3 



2 I 



jn + q) 



« 2 IS 



2i ( n - q) 



Sill" I — I SI IT | 7T 



(2.28) 



which can be directly obtained from l|2.15l) by replacing momenta with their lattice counter- 
parts k — > fc = - sin (f fc). 

The lattice quantum action is (dropping again the eigenvalue A3 since it is independent 
oiq): 



N t -1 r l 



s quL = y 

n=0 

+ log 



d 3 k' 



(2^) 3 



log 



siir 1 — I 1 snr I ?\ 



> + g) 
N 



,2 / ft i 



sm — + sm 7r 



Nt 



(2.29) 
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Figure 3: The lattice effective potential Vi(q, N t ) — Vj,(0, N t ) for different Nts, compared to the 
continuum effective potential q 2 (l — q) 2 . 

After performing the sum, the lattice quantum action is 

Sf lL = 2iV s 3 / log[l-2cos(27rg) e - 2A, *' l + e- 4A,f ' 1 ] (2.30) 

Jo l^ 7r ) 

h = log(K + VK 2 + 1) (2.31) 
K 2 = sin 2 (^) w(|)+sin 2 (|) (2.32) 
which is the same as stated by Weiss in \V2\ . 

We want to compare this to the continuum expression to make a quantitative prediction 
for corrections to the interface tension measured on the lattice. From (12.181) we have in the 
continuum (V — L x L y L z ) 

= lv"TV [q 2 (l~q) 2 ] (2.33) 



Defining 



V L (q,N t ) = [^N^j fj^E log[l-2cos(2^)e- 2 ^ + e - 4 ^] 



(2.34) 



and subtracting Vl(q = 0) from the action in order to make it zero when the system is in 
one of the vacua we get 

Sr L = \VT\ 2 (V L (q,N t )-V L (0,N t )) (2.35) 

Comparing this to we see that V(q) = q 2 {l - q) 2 is replaced by V L (q, N t ) - V L (0, N t ) 

on the lattice. Both are shown in Figure |21 

Taking the continuum limit Nt — > oo we have h — > ^ and therefore 

, i / 3 „,\ f°° 4nk 2 dk , r „ „ ,„ , 

VM = (^ 3 J I w -log[l-2oo 8 (2,rff)e 

= fr"?) y a^ 2 ^ l°g[l - 2cos(27rq)e 



-W t fe + e -2JV t fei 



-x ,_ e -2xl 
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N t 


Clat 


1 


1.2237 


2 


1.4233 


3 


1.29467 


4 


1.17025 


5 


1.09927 


6 


1.06278 


7 


1.04324 


8 


1.03181 


9 


1.02451 


10 


1.01953 



Table 1: Lattice correction factors Ci a t versus the number Nt of temporal slices. 



0.1 



0.01 



0.001 




relative correction 
1 .86 / Nt 2 



0.01 



0.1 



1/Nf 



Figure 4: The relative correction to the interface tension a, (Ci a t(Nt) — 1), for different values of 
Nt- Only for large values of Nt does the relative correction show the expected behaviour oc 



which can be shown to be identical to q 2 (l — q) 2 
to q. 



constant, by differentiating with respect 



This gives for the lattice correction factors Ci a t(Nt) = a(N t )/a(N t — oo): 



Q at (N t ) = 6 I dq^V L (q,N t )-V L (0,N t ) 
o 



(2.36) 



A numerical evaluation yields the factors in Tabled As is to be expected in a bosonic theory 

the corrections are of order 0(a 2 ), i.e. the factors are of the form ^1 + 0(^)j. This is 

shown in Figure 0] However, this asymptotic behaviour sets in only for large values of Nt. 
For small values of Nt, the ones accessible to numerical simulations, the correction factors 
show a quite different, not even monotonic (N t = 1), behaviour. This knowledge is crucial 
when trying to extrapolate to the continuum limit with only small- N t Monte Carlo data. 



To leading order the lattice correction factors C'i a t do not depend on the number of col- 
ors N, which only occurs as a prefactor in the 1-loop quantum action and cancels out when 
comparing the continuum to the lattice effective potential. Thus Table I applies to SU (N) 
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as well. 



There is one more thing to say, concerning finite size effects. So far, we chose the spa- 
tial dimensions N s to be oo. But we can also calculate the impact of a finite size in the 
spatial dimensions on the lattice correction factors. Let us consider a cubic geometry for 
simplicity. For finite N s the integral over the momenta k becomes a sum and we obtain 

„ //V\ 3 N '- x 

V L (q,N t ,N s ) = £ log[l-2cos(2^)e- 2 ^ + e- 4 ^] 

K 2 = sin 2 (jj-n^j + sin 2 f^-%J + sin 2 fe»»J ( 2 - 37 ) 

One finds for example at N t — 2 and aspect ratio p = 2 (p = j£ ) that Ci a ± in Table 
should be multiplied by 0.885, for p = 4: 0.966 and for ,0 = 8: 0.992*. This shows, that with 
increasing aspect ratio, the impact of a finite N s becomes quickly negligible. If finite size 
effects are taken into account, Ci a t decreases slightly, and therefore the corrected interface 
tension a(T)/Ci a t increases by a tiny amount. 



3 Comparing with Simulations 

3.1 Measuring the Interface Tension 

A direct Monte Carlo measurement of a from the free energy definition (|1.4I1.2|) is imprac- 
tical. First, there is an overlap problem: configurations contributing to the numerator and 
the denominator of (|1.4|l are physically, macroscopically different, since they respectively 
contain an interface or don't. Thus, importance sampling with respect to the denominator 
fails for all but the smallest volumes. Second, the interface is translationally invariant, so 
that the corresponding entropy must be carefully subtracted to avoid logL corrections. Both 
problems can be solved by a clever choice of algorithm. Historically, this choice has been 
arrived at in two steps. 

An elegant approach to the overlap problem is provided by the 'snake' algorithm |SJ|2|. 
It consists of building the interface in N x x N y steps, one plaquette at a time, and factorizing 
the ratio ea. (jl.4j> : 

Ztbc _ Zjy Zn-1 Z\ ^ 

Zpbc Zpj-i Zn-2 Zq 

where N — N x x N y , and Zk, k = 0,...,N (Zp? = Z t b c and Z$ = Z p b c ) is the partition 
function for a 't Hooft loop of area ka 2 , with a stack Vk of k corresponding plaquettes 
"twisted", i.e. having the coupling (3 changed to —(3: 

Z k = [ PC/cxp(+/3E( 1 - T^RcTr U P ) - pj^i 1 - ^ cTr U P )) (3.2) 

where Vk is the complement of Vk- When mod (fc, N y ) — 0, the 't Hooft loop is rectangular 
with perimeter 2N y , otherwise it has perimeter (2^ + 2), and one side contains 2 kinks (the 
bold x-links in Fig. [SJ . The progressive flipping of the plaquettes is described by a snake- like 
motion, hence the name attached to the algorithm 9 . Each ratio in l|3.1|) is of 0(1) and can 
be efficiently estimated by a separate Monte Carlo simulation. 

In fact, this algorithm can be simplified 10 , and the N x x N y Monte Carlo simulations 
reduced to a single one, by the following observation. — ^ s a ra ti° °f t w0 't 

Hooft loops differing by a 2 in area. If mod (k, N y ) ^ {0, —1}, they have the same perimeter 
and number of corners. They form the dual analogue of a Creutz ratio, commonly used to 
measure the ordinary string tension. Therefore, 

lim — log fc+1 = aa 2 . (3.3) 

JVx.JVy^OO Zk 
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Figure 5: Illustration of the 'snake' algorithm, which builds the interface one plaquette at a 
time, and of its simplified version, which requires a single Monte Carlo simulation of a system 
with a partial interface. 

A crucial advantage over the Creutz ratio is that the integrand in Zk is of exponential form: 
it is positive and can be used as a sampling probability. Then, a given statistical precision 
on Zk+i/Zk requires the same computer effort, independently of the 't Hooft loop size. As 
stressed in |13) . this is in sharp contrast to the ordinary string tension, where the statistical 
accuracy degrades exponentially with the Wilson loop area. Therefore, (|3.3|1 can be rewritten 
and measured as an expectation value with respect to Z^: 

aa 2 = Urn -log (exp^/^ReTx U Pk+1 )) (3.4) 

-Vc ,iV y — >oc Z 

In our simulations, we use the equivalent (exp(— /J^ReTr [/p fe ))o/(exp(+/3|ReTr Up k+1 )) , 
where the expectation value (..)o is taken in the ensemble which interpolates between Zk 
and Zk+i, where the coupling attached to Up k+1 is set to zero. This observable has a smaller 
variance, which can be further reduced by performing multi-hit (analytically) over the 4 links 
making up Up k+1 , since they are now decoupled from each other. 

Note that translation invariance is broken. Our 't Hooft loops Zk and Zk+i have fixed 
boundaries, and the observable is a function of a single plaquette in the system. This solves 
the second problem mentioned at the beginning: there is no entropy to subtract from the 
free energy of a partial interface. 2 It also suggests to organize the Monte Carlo updates 
hierarchically in shells centered around Pk+i ^3]. The gauge links of the innermost shells, 
upon which the observable depends most sensitively, are integrated more thoroughly with 
more frequent Monte Carlo updates. 

Since taking a large 't Hooft loop in (|3.4|) requires no more work than taking a small one, 
we choose the size as large as practical, as per Fig. El with length N y and width r ~ N x /2. 
r can be increased even more, but taking it too close to N x brings about the possibility of 
reducing the total free energy by trading the single interface of width r for two of them: one 
interface of width [N x — r), and one full interface of width N x . Translation invariance of 
the latter causes the free energy reduction. The choice r ~ N x /2 completely eliminates this 
finite-size effect. 

2 In the originai 'snake' algorithm, translation invariance is restored in the last ratio Zn / Zn-i- Its Monte 
Carlo evaluation is plagued by very long autocorrelation times, corresponding to translations of the full interface. 



12 



For large 't Hooft loops of size N y , r £ = 1/ \Jaa 2 , the leading deviation of (|3.4|l from its 
asymptotic value comes from two contributions: Gaussian fluctuations of the interface away 
from its minimal area (known as capillary waves, or Luscher correction), and interaction 
between the two kinks in the 't Hooft loop perimeter Fig.^f. (i) The Luscher correction can 
be evaluated analytically. For a vibrating L y x r interface with £ <C r <C L y , it takes the 
familiar form a — > a + t^-tt- In our case where r ~ i-^j/, it is obtained from the partition 
function appropriate to our boundary conditions |14| 



7 — 2 

'-'surface — ^] 



2r 



t,(t) = e^ T JJ (1 - e i27rnr ) , (3.5) 



71=1 



We subtract this analytic correction from our Monte Carlo results, (ii) A kink of size a in the 
't Hooft loop perimeter creates UV excited string states of energy 0(a^ 1 ), which propagate 
along the y direction over a distance s rs aN y /2. Therefore, we expect a correction of 
order exp(— N y /2) due to the two kinks. Indeed, Monte Carlo data obtained on lattices of 
increasing size L 3 x 2 at fixed j3 are well described by the ansatz (aa 2 + Luscher correction 
+ci exp(— C2-/V y /2)) with ci,C2 ~ 1 ^H]. This finite-size correction is negligible compared to 
our statistical errors provided N y ^> — 21og(cra 2 ). Numerical tests at several values of (3 and 
several lattice extents Nt indicated that finite-size effects cause an underestimate of the true 
interface tension, and led us to aspect ratios {N x ,N y ,N z }/Nt as large as 16. 

To increase statistics, we typically estimated 2 ratios Zk+i/Zk, corresponding to plaque- 
ttes Up k+1 at the "center" of the lattice (i.e. with coordinates (^S + {0, 1})). For each 
run, an average of 20k multi-hit, multi-shell measurements were taken. 



3.2 Results versus 2-Loop Perturbation Theory 

In Section[3 we derived the interface tension at 1-loop, in the continuum and on the lattice. 
The continuum 2-loop result has been obtained in |S] and in 1 1 1 1 , where the running coupling 
g(T) is chosen such that the 1-loop effects disappear in the renormalization of the coupling 
in the dimensionally reduced theory. For SU(N): 



with Euler's constant = 0.577215.... The coupling g(T) is given in ^5] as 

g 2 (T) = 9 y(ri 

which defines the A-parameter 



M =477Te-(TE + i/22) 



(3.7) 



e ( 7 E+l/22) 
47T 



A = JZ A ms (3-8) 



With this result it is possible to express the interface tension in (|3.fi|l to subleading order in 
terms of the lattice bare coupling g^AriT). To do that, we use the SU(2) relation 

A-^ = 19.8228A Lj4 t (3.9) 

to obtain the 1-loop result 

9 2 = 9 2 lat^ + 0.10016 S 2 ) (3.10) 



3 Excited states of the interface also bring a correction. Assuming an energy gap ir/r between the ground- 
and first excited state propagating along the y direction, one can check that this correction is smaller than our 
statistical errors, for all lattice sizes and couplings considered here. 
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Figure 6: The dimensionless ratio of the interface tension over the temperature squared, 
Tj%(Nt,/3), versus (3, for different values of N t . The data at high (3 are described almost perfectly 
by the parameter- free formula ()3.12j) . after the lattice correction factors Ci a t{Nt) have been 
inserted. For Nt = 1, the perturbative curve describes the measured data almost perfectly over 
the whole range (3 > f3 c . The reason for this remarkably good agreement is, that the singularity 
in (|3.12|) occurs for (3 = 0.8587, which is fortuitously close to the critical value f3 c = 0.8730(2) 
PS! of N t = 1. 



which is substituted into i|3.6|) to get 

. 4^ 2 T 2 
a(T) 



3V6 9lat{T) 



(1 - 0.2U67g 2 LAT (T)) 



(3.11) 



But we want the lattice bare coupling to run with the scale a -1 instead of T = (a-/V t ) _1 
because what is used in a simulation is (3 — 2N / g\ AT (a~ l ). This can be easily translated 
using the /3-function i|3.16[l . Also the lattice correction factors Ci a t(Nt), calculated in Section 
12.21 and given in Tabled need to be included. Altogether, one gets for SU(2) 



-&(l 
3V6 2 V 



(0.21467- 



■A)log(JVt)) ~ R 
P 



(3.12) 



to subleading order, where (3$ is given in l|3.15|l . Note that we have neglected 0(g 2 ) cor- 
rections in the lattice correction factors Ci at (N t ). How well this (parameter free!) formula 
describes the measurements of the interface tension at different A^'s for high (3 can be seen 
in Figure HO The N t = 1 comparison is shown separately because the data do not lie neatly 
above the other curves but cross right through. 



We have not performed a similar comparison for SU(3). But analogous steps predict: 

^(P,N t )=C lat (N t )^- (0-34805 + fi,log(JV t )) |) (3.13) 



3.3 Results versus Temperature 

Let us collect our measurements of the interface tension for different values of (3 and Nt , and 



correct them with the appropriate lattice correction factor Ci a t(Nt), forming -SziT) 



/corrected 



tjv(/3, N t )measured/Ci a t(Nt). The bulk of the cutoff effects has been removed with the division 
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Figure 7: The ratio of (f) corrected over ^(f )cKA (ea. (|3.17j) ) versus f = T/K-jjg. Deviations 
from 1 are small (~ 2%), and become appreciable only when T<WT C , with T c /A-^g = 1.31(8) 

CHI- 



by Ciat(Nt), so that the left-hand side can be compared with the corresponding continuum 
quantity evaluated in perturbation theory, as a function of temperature. Since we lack the 
possibility to set a physical temperature scale (to do that, we would need the non-perturbative 
/^-function obtained, e.g., from the step scaling function J7j), we calculate the temperature 
T = l/(a((3)N t ), associated with each value of using the perturbative 2-loop formula for 
the lattice spacing 



-/3i/(2/3 2 ) 

a ^ ALAT) = A^ T ^{-W ){f ] (3 ' 14) 



1 ( (3 \ f4[3 

U(p,l±LAT) = 

with (N = 2 here): 



which is justified as long as T ^> A-p-g. We want to express the temperature in terms of the 
parameter Aj^-g. One obtains 

19.8228 / /? \ /4/3 ^ _/3l/(2/3 ° 2) 



From equation l|3.6[) . we know the relation between ^ and the coupling g at 2-loop order. In 
|2], Giovannangeli and Korthals Altes show that higher order corrections to equation (|3.fci|) 
are very small. Thus one may hope, that the evaluation of (|3.6ll for SU(2): 

a 4ir 2 1 

t^ {t)gka = ^7EW) (1 " °- 16459 ^ (T) ) (3 - 17) 

is sufficient to describe ^ over a large range of temperatures. For g 2 (T) we use the solution 
to the renormalization group equation to 2-loop: 

g 2 (T) = = (3.18) 

ft log(T2/A2) + 09a/^b) loglog(T2/A 2 ) 

where A is related to A-p^- through equation l|3.8|) . The perturbative approximation ^ (T)gka 
l|3.17|l is compared to our measurements ^(T) correcte d in Figured which shows the ratio 



15 



<?{T) CO rrected/o'{T)GKA versus T / Ajjg. As can be seen, the deviations from 1 are remarkably 
small and remain at the 2% level down to temperatures O(10)T c . This figure is a beautiful 
confirmation of the work of Giovannangeli and Korthals Altcs. 

4 Discussion 

We have shown that cutoff effects on lattice measurements of the Yang-Mills interface tension 
(or dual string tension) are large for practical lattice sizes used in Monte Carlo simulations. 
They do not disappear as f3 — > oo as could be naively expected. Table I, which lists these 
corrections calculated at 1-loop order as a function of the temporal lattice size N t for any 
SU (N) group, should be of practical use. 

Moreover, we have shown precise agreement between our SU(2) Monte Carlo data and 

2- loop continuum perturbation theory, once the correction factors above have been folded in. 
The ratio of the measured over the perturbative interface tension remains 1 within ~ 2%, 
from very high temperatures down to O(10)T c . Our results put the perturbative calculation 
of the interface tension on firm ground, and show that the perturbative expansion converges 
well for this quantity, in contrast to the notorious expansion of the pressure . The reason 
for this remarkable contrast lies presumably in the topological nature of the 't Hooft loop. 
It naturally encodes the Zjq center symmetry, which is absent in the perturbative expansion 
of the pressure. 

The small remaining discrepancy between the measured data and perturbation theory 
can be assigned various causes. We list them in a subjective order of decreasing importance: 
(i) 0(g 3 ) correction in a 3 ; (ii) 0(g 2 ) corrections in the lattice correction factors Ci a t\ {Hi) 

3- loop terms in the lattice /3-function 1)3. f 4(1 . which will change the £-axis T/Aj^g in Figure 
[3 (iv) 3-loop terms in the running of gr(T) 1(3. f 8|l : (v) remaining finite size effects in the 
Monte Carlo simulations. 

Of course, the low-temperature drop in Figure \7\ is natural. The interface tension is 
a (dual) order parameter for the phase transition between the cold confining and the hot 
deconfining phases, which for SU (2) is second-order. The corresponding correlation length 
therefore diverges near the critical temperature T c like (T/T c — , with the critical expo- 
nent v of the 3d Ising model as expected from universality jHJ. This singularity is not seen 
by the perturbative calculation. Hence, the ratio (7(T) correcte d/&(T)cKA in Fig. goes to 
zero like {T/T c - l) 2v . 

It would be useful to extend our calculation of the lattice correction factors Ci a t to order 
g 2 , for arbitrary SU{N) group. At this order, one may observe a dependence on the relative 
Zn separation of the two vacua on either side of the interface. This dependence would enter 
in the comparison of dual fc-string tensions (ak/^i){T) between lattice measurements |1(J) 
and continuum perturbation theory 

Finally, we became aware of Ref.JUj as this work was being written up. Ref. deals 
with the same issue of connecting lattice measurements and continuum calculations of dual 
tensions, for the 3d theory (see Table VIII there). Nevertheless, it seems that some of their 
analytic results generalize to the Ad case we studied. 

Note added 

The generalization mentioned above has just been presented in Ref. [23]. Table I there 
contains similar lattice correction factors as our Table I. The small numerical differences can 
presumably be assigned to truncation errors in the series used in Ref. |2U| . The numerical 
simulations in the paper focus on ratios of fc-interface tensions, or rather of their derivatives, 
in SU(N), N > 3. 
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